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ABSTRACT 

We present an analysis of X-ray and UV data obtained with the XMM-Newton 
Observatory of the long period dwarf nova RU Peg. RU Peg contains a massive white 
dwarf, possibly the hottest white dwarf in a dwarf nova, it has a low inclination, thus op- 
timally exposing its X-ray emitting boundary layer, and has an excellent trigonometric 
parallax distance. 

We modeled the X-ray data using XSPEC assuming a multi-temperature plasma 
emission model built from the MEKAL code (i.e., CEVMKL). We obtained a maxi- 
mum temperature of 31.7 keV, based on the EPIC MOS1, 2 and pn data, indicating 
that RU Peg has an X-ray spectrum harder than most dwarf novae, except U Gem. 
This result is consistent with and indirectly confirms the large mass of the white dwarf 
in RU Peg. The X-ray luminosity we computed corresponds to a boundary layer lu- 
minosity for a mass accretion rate of 2 x 10~ 11 Mq/yt (assuming M W( j = 1.3Mq), in 
agreement with the expected quiescent accretion rate. The modeling of the O viii emis- 
sion line at 19A as observed by the RGS implies a projected stellar rotational velocity 
v ro ts'mi =695 km s _1 , i.e. the line is emitted from material rotating at ~936-1245 
km s _1 (i ~ 34° — 48°) or about 1/6 of the Keplerian speed; this velocity is much 
larger than the rotation speed of the white dwarf inferred from the Far Ultraviolet 
Spectroscopic Explorer (FUSE) spectrum. Cross-correletion analysis yielded an unde- 
layed (time lag ~ 0) component and a delayed component of 116±17 sec where the 
X-ray variations/fluctuations lagged the UV variations. This indicates that the UV 
fluctuations in the inner disk are propagated into the X-ray emitting region in about 
116 sec. The undelayed component may be related to irradiation effects. 

Subject headings: accretion, accretion disks, binaries: close, Stars: white Dwarfs, Stars: 
dwarf novae (RU Peg), X-rays: binaries 

1. Introduction 

1.1. The long period dwarf nova RU Peg 

Dwarf novae (DNe) are a class of weakly-magnetic cataclysmic variables (CVs) which are 
interacting compact binaries in which a white dwarf (WD, the primary star) accretes matter and 
angular momentum from a main (or post-main) sequence star (the secondary) filling its Roche-lobe. 
The matter is transferred, at a continuous or sporadic accretion rate (M), by means of an accretion 
disk usually reaching all the way to the WD surface. Ongoing accretion at a low rate (quiescence) 
is interrupted every few weeks to months by intense accretion (outburst) of days to weeks (a dwarf 
nova accretion event). DNe are powerful X-ray sources with luminosities of 10 30 — 10 33 erg s -1 . 
The X-ray emission is believed to originate in the boundary layer (BL) between the slowly rotating 
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accreting WD and the fast rotating (Keplerian) inner edge of the accretion disk, where the material 
dissipates its remaining rotational kinetic energy before it accretes onto the surface of the WD. The 
typical DNe, i.e. those systems ex hibiting normal DN ou tbursts are the U Gem sub-type of DNe 
(according to the classification of iRitter fc Kolb I (|2003l ) ) and are located above the period gap 
(P orb > 3hr). 



RU Peg is a U Gem type DN with an orbital period P or b 
K2-5V, a primary (WD) mass M„, rf = 1.29+ 016 



O 2o^0) and a secondary mass Mi 



8.99 hrs, a secondary spectral type 
0.94 ± O.O4M 



(jStover lll98ll ; IWade lll982l ; IShafter Ill983l ). The system has a magnitude range V ma x — Vmin ~ 9.0- 
13.1 with outbursts lasting ~20 days and recurring every ~50 days. The near-Chan drasekhar WD 



mass has been corroborated by the sodium (8190A) doublet radial velocity study of lFriend et al 



(|l990l ). They obtained a mass of 1.38±0.06 Mq for the WD and also found a range of inclination 
angles betwee n 34°-48° in agreement with the range of plausible inclinations found in the study by 
Stover I (|198ll ). More recently, a Hubbl e Fine Guidance Sensor (FGS) parallax of 3.55 ± 0.26 mas 
was measured by iJohnson et al. I (J2003J) yielding a distance of 282 ± 20 pc. 



RU Peg was observed with IUE under several different obser ying programs both in 



quiescence 



and in outburst and was par t of several survey-like studies (e.g. lLa Dous et al. I (|l985l ); IVerbunt 
(|1987l ): ISzkody et al. I ()199ll ) to cite just a few). 



Sion fc Urban I (120021 ) modeled 4 IUE spectra obtained in deep quiescence with accretion disks 
and photospheres. They found that a very hot WD dominated the FUV spectrum with a temper- 
ature T e ff = 50 — 53, 000K which places RU Peg among the hottest WDs in DNe. The distance 
corresponding to their best fitting, high gravity (Log(g) = 8.7) photosphere models was 250 pc. 



More recently, (IGodon et al. 1 120081 ) modeled the FUSE spectrum of RU Peg in quiescence 



and obtained a WD with a temperature of 70,000K, a rotational velocity of 40 km s -1 , assuming 
Log(g) = 8.8 and a distance of 282 pc. In this later study the higher temperature obtained in the 
model fitting is mainly a consequence of the assumed larger distance and gravity. It is clear that 
RU Peg has a massive WD and therefore a deep potential well, and its surface temperature is very 
large (>50,000K), possibly pointing to strong accretional and boundary layer heating. For these 
reasons, we chose RU Peg as our X-ray target, as it is expected to be a copious source of X-rays 
and should be an ideal candidate to study its boundary layer. 



1.2. The Boundary Layer 

The boundary layer is that region between the slowly rotating accreting WD and the fast 
rotating (Keplerian) inner edge of the accretion disk. For accretion to occur, gravity has to overcome 
the centrifugal acceleration, and this happens in the boundary layer as the material dissipates its 
remaining rotational kinetic energy before being accreted onto the surface of the WD. In the 
boundary layer the rotational velocity is sub-Keplerian and decreases inwards. 
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to the disk luminosity (jPringle 



The standard disk theory rtShakura Sunvaev 1119731 ) predicts a BL luminosity almost equal 



1981 



l bl = (i - f ) L disk = (i - n 



2 GM*M 



2R* 



(1) 



where G is the gravitational constant, M* is the mass of the accreting star, R* its radius, M is 
the mass accretion rate, and f3 is the stellar angular velocity in Keplerian units f3 = f2*/0#(i2*). 
For most systems the WD rotational velocity is of the order of a few 100 km s _1 and one therefore 
has j3 « 1 and Lbl ~ Ldisk- For a star rotating at (e.g.) 10% of the breakup velocity eq.(l) 
gives Lbl = 0.99L,fefc- However, the standard disk t heory does not ta ke into account the boundary 
layer explicitly. Using a one-dimensional approach, iKluzniak I (|1987l ) has shown that part of the 
BL (kinetic) energy actually goes into the spinning of the star in the equatorial region, and the BL 
luminosity is more accurately given by the relation 



Lbl = (1 — /3) 2 Ldisk- 

For a star rotating at 10% of the breakup velocity, eq.(2) gives Lbl 
cantly different than eq.(l). 



(2) 

0.81Ldiski which is signifi- 



Because of its small radial extent, the BL is expected to emit its energy in the X-ray bands 
(Lbl ~ L x-ray)- At high accretion rate (or durin g dwarf nova outburst), the BL is expected to b e 



optically thick with a temperature T ~ 10 5 — 10 6 K ( Godon et al. Hl995l ; |Popham fe Naravan 



and emits in the sof t X-ray band. Observations of CVs in high state (e.g. iMauche et al. I ([1995); 



1995) 



Baskill et al. I (J2005J)) confirm these p redictions. Accretion onto white dwarfs in some symbiotics 



prese nts similar soft X-ray emissions (|Luna &: Sokoloski 1 120071 ; iLuna et al. 1 120081 ; iKennea et al 



2009) because of the large mass accretion rates in these systems. At low mass accretion rate (or 
during dwarf nova quiescence), as the density is muc h decreased, the BL becom es optically thin 
and emits in the hard X-ray band with T ~ 10 8 K ( Naravan fc Popham I Il993l ) (see below for 
observational evidence of such optically thin BLs). Hence, during quiescence the emission should 
arise from a very hot plasma very close to the WD surface. 



Pr evious X-ray observations of dwar f nova systems in quies cence (e.g. Ivan der Woerd &: Heise 



(|1987l ); Ivan Teeseling &: Verbunt I (|l996l ); iBelloni et al. I (|l99ll )) while confirming the presence of 



hard X-ray, deduced that, contrary to the theory, the quiescent BL luminosity was underluminous. 
Namely, they found Lbl = Lx-ray « Ldisk, which confirmed the original claim that BLs were 
actually missing (Ferland et al. 1982). However, these earlier results assumed the disk to be 
the source of the optical and ultraviolet radiation (Ldisk = L opt + Lrjy) and used eq.(l) rather 
than eq.(2). Mo re recently, X- r ay Multi M irror-Newton (XMM-Newton) observations of 8 DNe 
in quiescence by iPandel et al. I ([2003, 1200a ) revealed that a significant part of the emitted FUV 
flux (Ljjv) actually originates from the WD itself, and the e vidence for underlum inous BLs in 



quiescent DNe was refuted. For RU Peg, FUSE observations (IGodon et al. 1120081 ) indicate that 



the WD contributes possibly most of the FUV light with a temperature T > 50, 000K. In addition, 
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Godon Sion I (120051) noted th at the region where the boundary layer meets the inner edge of the 
Keplerian disk (jPopham II 19991 ) can also contribute some FUV flux. It is clear now that one cannot 
just compare the X-ray luminosity to the optical + UV luminosity to check whether the boundary 
luminosity is as large as the disk luminosity. 



For the 8 dwarf novae caught in quiescence, iPandel et al 



(j2003l . 12005J), using the XMM- 
Newton data, obtained X-ray boundary layer luminosities of the order of ~ 1 x 10 31 ergs s _1 to 
6.6 x 10 32 ergs s _1 , with temperatures ranging from ~ 8 to 55 keV, and mass accretion rates deduced 
from X-rays in the range W~ 12 M@/yic to lO~ lo M0/yr. In this work, we present XMM-Newton 
observations of RU Peg taken in quiescence to derive its X-ray luminosity and gain information on 
its boundary layer. 



2. Observations and Analysis 



The XMM-Newton Observatory (iJansen et al. 1 120011 ) has three 1500 cm 2 X-ray telescopes 
each with an European Photon Imaging Ca mera (EPIC) at the focus. Two of the telescopes have 
Multi-Object Spectro meter (MOS) CCDs (jTurner et al. I hoOll ) and the last one uses pn CCDs 



Striider et al. 



2001 



(den Herder et al. 



2001 



for data recording. There are two Reflection Grating Spectrometers (RGS) 
The Optical Monitor (OM), a photon counting instrument, is a co- 



aligned 30-cm optical/UV telescope, providing for the first time the possibility to observ e simulta- 



neously in the X-ray and optical/UV regime from a single platform (jMason et al. 



2001 



)• 



RU Peg was observed (pointed observation) with XMM-Newton for an duration of 53.1 ks on 
June 9th, 2008, at 07:16:50.0 UTC (obsID 0551920101). At that time the system was at a visual 
magnitude of ~ 12.5, about 2 months into quiescence and 2 weeks before the next outburst (from 
AAVSO data0). Data were collected with the EPIC MOS and pn cameras in the prime partial 
window2 and prime full window imaging mode, respectively, the Reflection Grating Spectrometer 
and the Optical Monitor using the fast imaging mode (<0.5 sec time resolution) with the UVW1 
filter (240-340 nm). 

We analysed the pipeline-processed data using Science Analysis Software (SAS) version 9.0.0. 
Data (single- and double-pixel events, i.e., patterns 0-4 with Flag=0 option for pn and patterns < 
12 with Flag=0 for MOS 1,2) were extracted from a circular region of radius 60" for pn, and 45" 
MOS1 & MOS2 in order to perform spectral analysis together with the background events extracted 
from a source free zone normalized to the source extraction area. We checked the pipeline-processed 
event file for any existing flaring episodes and no sporadic events in the background were detected 
with count rate higher than 0.08 c s _1 (for MOSl,2) and 0.5 for pn detectors. Table 1 displays the 
background subtracted count rates for the EPIC pn, MOS1 and MOS2. 



http://www.aavso.org/ 
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The RGS observations were carried out using the standard spectroscopy mode for readout. We 
reprocessed the data using the XMM-SAS routine RGSPROC. We first made event files and deter- 
mined times of low background from the count rate on CCD 9 (which is closest to the optical axis). 
The final exposure times and net count rates showed that there were no sporadic high background 
events in our data. Table 1 displays the background subtracted count rates for RGS1 and RGS2. 
Source and background counts for the RGS were extracted using the standard spatial and energy 
filters for the source position, which defines the spatial extraction regions as well as the wavelength 
zero point. 



2.1. X-ray and UV Light curves 

The UV and EPIC pn X-ray light curves are shown together for comparison in Figure 1, 
where we have scaled the UV count rate to fit the X-ray count rate for a better comparison. 
The count rate for the UV data ranges between ~50 c/s and ~100 c/s with a time average of 
~71.8 c/s. For UV-bright objects, a count rate of 1 s _1 in the UVW1 filter translates into a flux 
of 4.5 x lCP 16 ergs cm - 2 s -1 A — 1 at 290 nm (see the online XMM-Newton documentation,.) . This 
gives for RU Peg a flux of ~ 2.2 — 4.5 x lCP 14 ergs cm - 2 s~ 1 A — ^ corresponding to a luminosity of 
r* 2.1 — 4.3 x 10 29 ergs s _1 A _1 (at a distance d=282pc). The time modulation of the UV data 
follows closely the time modulation of the X-ray data except around t ~12 ks, 20-22 ks, 31 ks, Sz 
48.5 ks, where the UV has relatively more flux for a duration of several hundred seconds (and up 
to 1,000 s). Since the UV is expected to be emitted further out than the X-rays, these four epochs 
where the UV light curves do not decrease as much as the X-rays might be due to the occultation 
of the X-ray emitting material by the WD, while the UV emitting region is not hidden from the 
observer. 

In order to study the correlation between the X-ray and the UV variability, we calculated the 
cross-correlation between the two light curves. We used time bins of 5 sec averaging several power 
spectra with 128 bins for the analysis. The resulting correlation coefficient as a function of time lag 
is shown in Figures 2a &; 2b. The correlation coefficient is normalized to have a maximum value 
of 1. The curve shows a clear asymmetry indicating existence of time delays. We, also, detect a 
strong peak near zero time lag suggesting a significant correlation between X-rays and the UV light 
curves. We call this the undelayed component (see Figures 2a & 2b). The positive time-lag in the 
asymmetric profile shows that the X-ray variations are delayed relative to those in the UV. In order 
to calculate an average time-lag that would produce the asymmetric profile, we fitted the varying 
cross-correlation by two lorentzians, with time parameter fixed at 0.0 lag and the other set as free. 
The resulting fit yields a lag of 116±17 sec. This is the delayed component. 



2 



http://xmm.vilspa.esa.es/external/xmm_user_support/documentation/uhb/index.html 
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2.2. EPIC Spectrum 

We performed spectral analysis of the EPIC data using the SAS task ESPECGET and derived 
the spectra of the source and the background together with the appropriate response matrices and 
ancillary files. How the photons were extracted is described in Section 2. The EPIC pn, EPIC 
MOS1 and EPIC MOS2 spectra were simultaneously fitted to derive the spectral parameters. The 
spectral analysis was performed using XSPEC version 12.6.0q (Arnaud 1996). A constant factor 
was included in the spectral fitting to allow for a normalization uncertainty between the EPIC 
pn and EPIC MOS instruments. We grouped the pn and MOS spectral energy channels so that 
there is a minimum of 80 (MOS1,2)-150 (pn) counts in a bin to improve the statistical quality of 
the spectra. The fits were conducted in the 0.2-10.0 keV range. The simultaneously fitted spectra 
from the three EPIC instrum ents are shown in F igure 3. We modeled the X-ray spectrum of 



RU Peg in a similar fashion as IPandel et al. I (|2005l ) and fitted the data with (TBabsxCEVMKL) 



model within XSPEC. TBabs is the Tuebingen-Boulder ISM absorption model (Wilms, Allen and 
McCr ay 2000) and CEVM KL is a multi-temperature plasma emission model built from the mekal 



code (IMewe et al. Ill985l ). Emission measures follow a power-law in temperature (i.e. emission 
measure from temperature T is proportional to (T/T max ) a-1 ). The residuals in Figure 3 show 
systematic fluctuation around the 6.7-6.9 keV iron line complex mainly from the EPIC pn data and 
some small low energy fluctuations exist in the MOS2 data, as well. This is due to the CTI (Charge 
transfer inefficiency) problem, in the pn (and possibly MOS2) instrument. These generally occur 
around lines due to small calibration errors and mostly effect only the line shapes leaving systematic 
residuals and increasing the reduced x 2 °f the fits. Our EPIC MOS1 data does not exhibit any 
CTI effects and the reduced \ 2 for the fit to this data alone is 1.15 (d.o.f. 290). The reduced x 2 
for the simultaneously fitted spectra is higher than the value for the MOS1 fit, but the spectral 
parameters for all three instruments are almost the same within the errors. Table 2 contains the 
spectral parameters from fits (using three detectors simultaneously) with the (TBabsxCEVMKL) 
model. Errors are given at 90% conf. level. We find a maximum plasma temperature in a range 
29-33 keV and mostly solar abundances of elements aside from oxygen and neon which we calculate 
to be subsolar. The unabsorbed X-ray flux is ^lj^xlO -11 erg s _1 cm -2 which translates to a 
luminosity of 4.1^Q3Xl0 32 erg s _1 at 282 pc (see sec. 1). The neutral hydrogen column density is 
4.3-4.5xl0 20 cm" 2 . 



2.3. RGS Spectrum 

The RGS analysis tool RGSPROC was used to obtain RGS1 and RGS2 spectra and to produce 
fluxed spectrum (i.e., using rgsfluxer) The resultant fluxed spectrum of the combined RGS1 and 
RGS2 detectors is shown in Figure 4 with line identifications. The detected lines and corresponding 
wavelengths are listed in Table 3. Fluxed spectra are obtained just by dividing the count spectrum 
by the RGS effective area. It neglects the redistribution of monochromatic response into the 
dispersion channels. Since proper response is not utilized by the fluxed spectrum, in order to 
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perform a similar spectral fit/analysis using the same (TBabsxCEVMKL) model within XSPEC, 
we used the count rate spectra produced for RGS1 and RGS2 simultaneously with the appropriate 
response files for each detector. This is a very efficient approach to find the spectral parameters 
in comparison with the EPIC results. The fitted RGS1 and RGS2 spectra are shown in Figure 5. 
Table 4 contains the spectral parameters from the fit with the (TBabsxCEVMKL) model within 
XSPEC. Errors are 90% conf. level and the fit is performed between 0.2-2.5 keV. 

The maximum temperature from the fits with the RGS data yield 24 keV with a large error 
range of 17-41 keV since the spectrum has a lower count rate and a lower upper energy boundary 
(i.e., 2.5 keV) compared with the EPIC data (i.e., 12 keV). We find that most of the RGS spectral 
parameters are consistent for the EPIC results. Oxygen and neon are subsolar in abundance 
and silicon appears to be slightly enhanced compared to solar abundance. The resolution of the 
RGS spectra is sufficient to measure the rotational velocity of the bounday layer via Doppler 
broadening of emission lines. We calculated the broadening using the O viii emission line at 19A 
which is the strongest line in the spectrum. We used a Gaussian model to calculate the a of the 
line (<t x 2.4=FWHM) along with a power-law/bremmstrahlung for the continuum. We find that 
FWHM=0.044A and the corresponding velocity in the line of sight is 695 km s _1 (AA/A = v/c). 
We also checked the resolution of RGS at around 19A and found 0.04-0.05 A and we caution that 
our measurement is on the limit of the RGS spectral resolution. 



3. Discussion 



The maximum X-ray shock temperature we obtain is 29-33 keV, based on the EPIC MOS1, 2 
and pn data, as they have higher count rates and braoder energy ranges than the RGS data. RU 
Peg has a harder X-ray spectrum than in most dwarf novae, but softer than U Gem i n quiescence, 



which also contains a massive white dwarf (ISion et al. 1 1 19981 : lLong fc Gilliland 1 11999). This result 



is consistent with and confirms the large mass of the white dwarf in RU Peg. 

The X-ray luminosity of RU Peg is 4.1 x 10 32 ergs s _1 (0.2-10.0 keV), and assuming we see 
only half of the boundary layer (if it is close to the star), we have Lu = 8.2 x 10 32 ergs s" 1 . Th is 



is in the range of the other quiescent DNe observed with XMM-Newton by lPandel et al. I (|2005l ) 



However, in order to fully compared RU Peg with these other systems, we also need to consider the 
UV luminosity. For the UV we use the spectral luminosity at 290 nm Ljjyw^ o btained from the 
OM data (see section 2.1). We have reproduced Figure 4 of iPandel et al. I (120051 ) in Figure 6 with 
the inclusion of RU Peg, which shows the quiescent DNe observed with XMM-Newton plotted on a 
Ljjvwi against L^i graph. In this figure the solid and dotted lines show the relation = L^sk for 
a UV luminosity predicted by a simple accretion disk model with an inner radius of 5000km and 
10,000km respectively. This simple disk model does not include UV contribution from the WD, or 
outer edge of the BL where it meets the inner disk. The location of a system in the vicinity of this 
diagonal (e.g. SU UMa, WW Hyi) indicates that the boundary layer luminosity Lm is comparable 
to the disk luminosity (here L^isk = Ljjvwi)- F° r PU Peg, as for U Gem for example, the excess 
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of UV is due to the contribution from the WD marked on the left of the graph. RU Peg has a 
UV luminosity corresponding to 53,350 K for a 8,000 km radius WD (assuming the flux scales 
simply as oc T 4 ), or 75,450 K WD for a 4,000 km radius (more consistent with the lar ge mass of RU 



Peg) . This confirms the temperature of the WD as derived from the FUSE spectrum (jGodon et al 



2008), and puts RU Peg (literally) in line with all the other quiescent DNe such as VW Hyi, U 
Gem, SU UMa, OY Car, and AB Dra, on the Luvwi versus Lu graph. 

Because of its large mass (and therefore small radius), the WD in RU Peg has a deeper 
potential well, and VKep of the order of 5-6,000 km s _1 in the BL. As the matter is decelerated in 
the BL, the X-rays emitting gas has velocities of a few 1000 km s _1 . The modeling of the O VIII 
emission line at 19A implies a projected rotational broadening v ro tsmi = 695 km s _1 , i.e. the line 
is emitted from material rotating at ~ 936-1245 km/ s _1 (since i ~ 34° — 48°) or about 1/6 of 
the Keplerian speed. This velocity is still much larger than the rota tion speed of the white dwarf 
inferred from the FUSE spectrum (40 km s^ 1 ; ( Godon et al. 1120081 )), This implies that the X-ray 
emission comes directly from the decelerating boundary layer material. It is possible that th e X-ra y 
emission originates in the equatorial region of the white dwarf as shown by lPiro &: Bildsten I (J2004J), 
who found that poloidal motion may be negligible at low mass accretion rates as a characteristics 
of DNe in quiescence. In such a case most of the dissipated energy i s radiated back into t h e disk , 
which may just ify the use of a one-dimensional treatment (such as iNaravan &: Popham I ([1993); 
Popham I (|l999l )). 



The X-ray luminosity we computed corresponds to a boundary layer luminosity (eq.2) for 
a mass accretion rate of 2.0 x 10~ n MQ/yr assuming a 1.29M Q WD mass, and it increases to 
3.0 x 10 _11 MQ/yr for a 1.2Mq WD mass. This is entirely consistent with a quiescent accretion 
rate. 

We checked the correlation between the variability of the X-ray data and the UV data and 
found two components; one delayed and the other undelayed. The X-ray and UV emission originate 
from distinct regions in the binary. Since RU Peg is a non-magnetic system, the X-ray emission 
is entirely from the inner edge of the boundary layer very close to the WD surface. As to the UV 
radiation, it is emitted by the heated WD, the very inner disk, and also that region where the 
outer boundary layer meets the disk (jPopham Ill999l ). The significant modulation correlation at 
Ai^O lag is expected to be caused by reprocessing of X-rays (i.e., irradiation by X-rays) in the 
accretion disk. Such time lags are on the order of milliseconds and proportional to light travel time 
which is well beyond the time resolution in our light curves. The delayed component of Ai~116 
sec lag is much longer and can not be produced by light travel effects nor by reprocessing of the 
X-ray, since the X-ray trails behind the UV. The only viable explanation, is that the time lag is the 
time it takes for matter to move inwards from the very inner disk (emitting in the UV) onto the 
stellar surface (emitting in the X-ray). This is the time it takes to spin down the material in the 
boundary layer T sp i n = At = 116s . The modulations of the UV component and its lagging X-ray 
counterpart are due to modula tions in M. A com parable time delayed (At~100 sec) component 
was also detected for VW Hyi (jPandel et al. 1120031 ) . and a much shorter one (~ 7s; Revnivtsev et 
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al. 2011) was detected for the intermediate polar EX Hya (indicating that the transit of matter 
through magnetic field lines to the poles is faster than though the non-magnetic boundary layer). 



Following the work of iGodon Sion I (|2005l ). we use the spin down time r sp i n = 116s and 
the rotation (or dynamical) time T ro t = 2nr /vk(i") =25. 3s to derive the viscous time in the BL 
t v = T^ pin /T rot =532s. This gives a boundary layer viscosity v = S>1i/t v = 4.8 x 10 13 cm 2 s _1 , 
where we have assumed M w d = 1.29M , R w d = 4, 000km , and a bounda ry layer size 8u given by 
the boundary layer radius ru = (1 + Sbi)R w d = 1-4-R^d (jPopham Ill999l ). In the alpha viscosity 
prescription v = ac s H, the value of the alpha parameter is then simply a = v/{c s H) ~ 0.003, 
where we ass umed c, 9 ~ 10 8 c m s _1 (for T ~ 10 8 K) and the vertical thickness of the boundary layer 



H « 0AR wc i (IPopham II1999T). This v alue o f a for the BL of RU Peg and the one derived for the BL 
of VW Hyi (0.009 - IGodon Sion I (|2005l )) indicate that the alpha viscosity parameter in the BL 
(otu) is much smaller than in the disk (a.dixk\ at least in quie scent dwarf novae) . T his is consistent 
with the analytical estimates of IShakura Sunvaev I (|1988l ) and IGodon I (119951) which obtained 
au << otdixk- Since the source o f the boundary layer viscosity is unknown (jlnogamow fc Sunvaev 



1999 



Popham &: Sunvaev Il200ll ). this result is important for future theoretical work investigating 



the source of the viscosity in the boundary layer. 



PG wishes to thank Bill Blair for his kind hospitality at the Johns Hopkins University where 
part of this work was carried out. Except for the symbol and temperature mark of RU Peg, 
Figure 6 was taken from lPandel et al. I (|2005l ). who kindly agreed that we reproduce Figure 4 from 
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Table 1: XMM-Newton Observations of RU Peg 



Instument 


Count Rate 




< > 


RGS2 


0.2881 ± 0.0030 


RGS1 


0.2351 ± 0.0029 


EPIC pn 


10.03 ±0.017 


EPIC MOS1 


3.154 ± 0.0084 


EPIC MOS2 


3.241 ± 0.0085 



Table 2: Spectral parameters of fit to the combined EPIC spectrum of RU Peg in the energy 
range 0.2-10 keV. Nh is the absorbing column, a is the index for power-law emissivity function 
{dEM = {T/T max ) a ' l dT/T max ) 

j T max is the maximum temperature for CEVMKL model. Element 
names stand for the abundance relative to solar abundances, Gaussian LineE is the line center for 
the emission line, uq is the line width; Kqevmkl and Kq are the normalizations for CEVMKL 
and Gaussian models respectively. The unabsorbed X-ray flux is given in the range 0.2-10.0 keV. 
All error ranges are given in 90% confidence level (Ax 2 =2.71 for a single parameter) 



PARAMETER VALUE 

N H (10 22 atoms/cm 2 ) 0.044+^1 

n 1 05 +0 - 03 

a i.uo_ 003 

T max (keV) 31.71^ 

O 3 +0 - 06 

Ne 0.55±g;^ 

Mg 1.3lij; 2 

Si 0.818:2 

S 0.9±g-| 

Ca 1.8±g : f 

Fe 0.8lgg 

Kcevmkl 0-0471q;qqJ 

Gaussian LineE (keV) 6.4 (fixed) 

o G (keV) 0.191^1 

K G 0.0000521°;™ 

Flux (10- 11 ergs c m - 2 S - 1 ) 4.ll°; 2 

Xl (d.o.f.) 1.5 (1741) 
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Table 3: Line Identifications of the XMM-Newton RGS spectrum of RU Peg 



Ion 


Wavelength 




< A > 


Mg xii 


7.20 




7.45 


Nex 


12.15 


Fe xvn 


15.03 


viii 


16.05 


Fe xvn 


16.80 


Fe xvn 


17.07 


viii 


19.00 


Ovn(i) 


21.80 


Ovn(f) 


22.11 


Nvn 


24.77 



Table 4: Spectral parameters of the fit to the RGS spectrum of RU Peg in the energy range 0.2-2.5 
keV with same definitions as in Table 2. 



PARAMETER VALUE 
N H (10 22 atoms/cm 2 ) 
a 

T max (keV) 

o 

Ne 
Mg 
Si 
S 

Ca 
Fe 

KcEVMKL 

Flux (10 -11 ergs cm~ 2 ) 
Xl (d.o.f.) 



031 + - 007 
u - uoi -0.007 

i 9 +0.16 

1 - z -0.16 

24 1+ 17 - 
7.0 

n e+o.13 
U.8_ 012 

i 1+0.5 
L - L -0.5 



2.T 



-1.2 
-1.1 
o 4 +2.0 



1.4* 



3.2 



1.0 (fixed) 
1.0 (fixed) 
u.u^o_ 005 



0.4 
0.4 



3.9 
1.3 (283) 
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Fig. 1. — The EPIC X-ray light curve (binned at 16 s, in red) together with the OM UV light curve 
(binned here at 100 s for display, in black; the data were binned at 20 s for the correlation check). 
The count rate of the UV data has been divided by 9.1736 to fit the count rate level of the X-ray 
data for easier comparison. The time t = corresponds to t = 3.29389320 x 10 8 s, counted from the 
MJD reference day 50,814, namely: 54,626.376 MJD. The time modulation of the UV data follows 
closely the time modulation of the X-ray data except around t ~12 ks, 20-22 ks, 31 ks, & 48.5 ks, 
where the UV has relatively more flux. 



-300' ' '-£00' ' '-\6ti ' ' 6 ' ' ' '1'00 ' ' '2'00 ' ' '300 
Time Delay (s) 



Fig. 2. — Cross-correlation between X-ray and UV light curves. Left: a single Lorentzian is used 
to fit the data at 0; Right: a second Lorentzian fit is added to fit the delayed component. The 
figures show strong correlation at zero time lag. The asymmetric profile centered at zero time lag 
(asymmetry towards positive lag) indicates that the X-ray variations are delayed with respect to 
the UV variations. Error bars indicate the standard deviations from the average of the value of the 
cross-correlation at any given lag in different time segments where power spectra are calculated. 
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Fig. 3. — EPIC M0S1,2 and pn spectra simultaneously fitted to the same model 
TBabsxCEVMKL. The model includes an ISM absorption model and a multi-temperature plasma 
emission model built from the MEKAL code. The maximum shock temperature obtained from the 
fit is 31.7 keV. The second panel below the fitted spectra is the residuals in standard deviations. 
One can notice a slight CTI problem of the pn detector around the iron lines (6.4-6.9 keV) where 
the fluctuations in the residuals due to CTI are evident. This does not affect the global fitting 
procedure but only slightly increases the reduced x 2 of the fits. 
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Fig. 5. — The simultaneously fitted RGS 1 and RGS 2 spectrum of RU Peg. The same model 
TBabsxCEVMKL is used for the fit. The residuals in the second panel below are in standard 
deviations. 
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Fig. 6. — The boundary l ayer (X-ray) luminosities vs. the UV spectral luminosity at 290 nm, as 



derived by iPandel et al. I (|2005l ) for quiescent DNe observed with XMM-Newton. We have added 
RU Peg on the graph. The solid and dotted lines show L&; = L^ s fc(= Luvwi) assuming UV 
luminosity predicted by a simple accretion disk model for an inner radius Ri n = 5000 km and 
10,000 km respectively. For systems higher above the line, such as RU Peg, the UV excess is due 
to the contribution from the WD. The UV luminosity of a WD with a 8000 km radius at various 
temperature is shown on the left. RU Peg corresponds here to a ~53,350 K WD with a radius of 
8000 km (which we have also added to the original graph), which translates to a ~75,000 K for 
a 4000 km WD, more in line with the large mass of the WD in RU Peg. WX Hyi was caught in 
outburst explaining its UV excess. 



